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Abstract 

I study a quantum Lifshitz point in a three-dimensional itinerant antifer- 
romagnet, in particular the scaling of the Neel temperature, the correlation 
length, the staggered susceptibility, the specific heat coefficient and the resis- 
tivity. At low temperatures, the model is shown to have the inverse staggered 
susceptibility and the resistivity varying as T 5//4 , and the specific heat coeffi- 
cient varying as T 1 / 4 . 

PACS numbers: 



1 



Typeset using REVTgX 



I. INTRODUCTION 



Magnetic metals close to a zero temperature phase transition have been a subject of 
active ongoing study§, mainly due to their apparent defiance of the Fermi liquid paradigm. 
Near a quantum critical point, one usually finds neither a nearly constant linear specific 
heat coefficient 7 = 5C /T ~ const., nor a virtually temperature-independent magnetic 
susceptibility x ~ const., nor a T 2 -resistivity p ~ const. + T 2 , all normally expected from 
a Landau Fermi liquid@. Instead experiments reveal a host of anomalous thermodynamic, 
magnetic and transport propertiesS calling for adequate theoretical description. 

In its current form, scaling theory of quantum criticality in itinerant magnets is due 
to Hert J and Millisi. It describes incipient magnetic ordering in an isotropic metal solely 
in terms of a boson order parameter with the correlation range in both space and time 
diverging at a quantum critical point. This divergence naturally leads to low-temperature 
anomalies in thermodynamic, magnetic and transport properties. Some of the early work 
on the subject was done by Mathoni, Makoshi and Moriyai, Ued J and Dzyaloshinskii and 
KondratenfeJi, who described magnetic transitions at T = using diagrammatic methods. 
These and other studies of nearly ferro- and antiferromagnetic metals as well as magnetic 
fluctuations in general, done in the spirit of self-consistent theory, were summarized in the 
book by Moriyai. A phenomenological approach aiming at establishing relations between 
various critical exponents at a quantum critical point has been developed by ContinentinoS. 

The Hertz-Millis theory and its implications describe rather well some itinerant ferro- 
magnets such as ZrZn2 and, to a lesser extent, MnSi. In ZrZn2, moderate external pressure 
p « 8 kbar is enough to reduce the Curie temperature Tq to zero. Near the critical pressure 
p c , T c scales as (p c -p) 3/i while the resistivity obeyffl 

p r*j j"i-6±o.i^ j n exce u en t agreement 

with the theoryiS'i. 

At the same time, several Ce-based antiferromagnets, where the Neel temperature can 
be tuned to zero by pressure or doping, appear to disregard the conclusions of the available 
theories. One of the most studied examples is0 CeCu 6 _, E Au r , where at the quantum critical 
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point the resistivity behaves as const. + T, and where over more than two decades in tem- 
perature the specific heat coefficient diverges as ln(T) while the inverse susceptibility fitsEl 
X' 1 ~ const. + T a with a » 0.8 ± 0.1. 

Another interesting case is CePd2Si2, a metal which orders antiferromagnetically below 
10 K. Applied pressure reduces the Neel temperature T/v from 10 K at the ambient pressure 
to about 0.4 K at 28 kbar. Between 15 kbar and 28 kbar, T/v falls linearly with pressure, and 
at 28 kbar CePd2Si2 becomes superconducting below 0.4 K, in a narrow strip between 23 
and 32 kbai§. At 28 kbar, the resistivity of CePd2Si2 exhibits a striking T 1 0-1 behavior 
over about two decades in temperature^. 

Neither CeCug^Au^ nor CePd 2 Si 2 appear to respect the available theoretical results 
which, for a bulk antiferromagnet, would imply the Neel temperature scaling! as T/v ~ 
(Pc — p) 2 ^ 3 , the specific heat coefficient!!! 7 ~ T 1 / 2 and the resistivity!'! p ~ T 3 / 2 . At the 
same time, theoretical results for two dimensions are much closer to what was observed in 
(fcOWta., as in this case the theory does yield const. + T resistivityS , logarithmically 
divergent linear specific heat coefficient and essentially linear scaling of the Neel temperature 
with pressure!. This led to an interpretation of data on CeCvLQ-. x Au x in terms of a purely 
two-dimensional magnetic ordering in an otherwise perfectly three-dimensional metal with 
only moderate anisotropySB. 

An alternative interpretation!^, based on the neutron scattering data, pointed at a pos- 
sible explanation in terms of highly anisotropic critical fluctuations with the stiffness in one 
of the directions vanishing at or very close to the quantum critical point. The conspiracy of 
this anisotropy with the anomalous frequency exponent^ a ~ 0.8 ±0.1 of critical fluctua- 
tions allowed to fit together the large body of experimental data obtained at the quantum 
critical point in CeCug^Auj;, including the specific heat, the uniform susceptibility and the 
neutron scattering scans. 

This latter scenario would imply a residual quartic dispersion (~ q 4 ) of critical fluctua- 
tions in the "soft" direction, thus also leading to the dimensionality reduction (by 1/2, as 
opposed to the purely two-dimensional order) and hence to qualitative modification of the 



Hertz-Millis theory. 

Finally, it was noticedQIll that the resistivity p ~ j"i-2±o.i Q f C e Pd2Si2, as well as p ~ 
j>i.25±o.i G f CeNi2Ge2 (which has the structure of CePd2Si2 with a smaller unit cell) may also 
be explained assuming anisotropic critical fluctuations with the residual quartic dispersion 
in one of the directions. 

All these experimental findings suggest that vanishing stiffness may be the aspect of 
physics needed to to describe the quantum criticality in Ce-based antiferromagnets by a 
Hertz-Millis type of theory. A critical point which, in addition to the onset of ordering, is 
characterized by disappearance of stiffness in one or several directions, is called a Lifshitz 
pointEl. In this paper, I study a quantum Lifshitz point, a curious yet possibly experimentally 
relevant coincidence of a quantum critical point (onset of ordering at T = 0) with a point 
where the stiffness vanishes in one or several directions in the momentum space. I develop 
scaling theory of a classical Gaussian regionB of a disordered phase near a quantum Lifshitz 
point in an itinerant three-dimensional magnet. I study a particular case of a Lifshitz point, 
where the incipient ordering exhibits anisotropic dispersion00§ which is quartic in only 
one direction and quadratic in the remaining two. 

The restriction to only the classical Gaussian region is due to the fact that in the other 
regions the low-temperature behavior is dominated by crossovers between various regimes 
and thus comparison to the experiment is rather hard to make. Moreover, the full phase 
diagram is sensitive to the relative strength of the coupling constants, whereas the results 
for the classical Gaussian region of a disordered phase are independent of this uncertainty 
and can be tested experimentally. 

It is important to note that Lifshitz point is a multicritical point and that, generally, 
one shall expect to find three phases in its vicinity, corresponding to one disordered and 
two different ordered states (see Section V). At the same time, none of the experiments 
on CeCue-xAuz or CePd2Si2 which I am aware of, indicated presence of more than one 
magnetically ordered phase. To that end, I show that, near a quantum Lifshitz point, one 
may indeed find only one ordered phase rather than two, which makes the model potentially 
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relevant to CeCvL 6 - x AuMM CePd 2 Si 2 and CeNi 2 Ge 2 §0§&i. 

Bearing in mind that both CePd 2 Si 2 and CeCu6_a;Au x are antiferromagnetic metals, one 
is led to assume that the critical mode is non-conserved and over-damped. Therefore this 
work may be viewed as an experiment-motivated extension of the Hertz-Millis theoryHS to 
a particular case of a quantum Lifshitz point. 

In Section II, I present the results and discuss them. I determine the equation of the 
critical line, the behavior of the order parameter susceptibility and the correlation length 
close to the quantum Lifshitz point. I also estimate the low-temperature behavior of the 
conductivity and the anomaly of the linear specific heat coefficient. Then I compare these 
theoretical results with the experimental findings@@0EHHl. 

In Section III, I derive the scaling equations, solve them in Section IV and obtain the 
results outlined in Section II. In Section V, I discuss the phase diagram near the quantum 
Lifshitz point and set the conditions for appearance of only one ordered phase. 

To demonstrate that the quantum Lifshitz point I study is well defined, in Section VI 
I examine how interactions generate stiffness in the "soft" direction, and show that this 
effect can be neglected at all relevant momenta and frequencies. Finally, in Section VII, I 
summarize the results and comment on them. The appendix provides the calculation details. 

Self-consistent treatment of a model with the stiffness vanishing in one of the directions 
in the momentum space has been given recently by C. Lacroix et a/.il with reference to 
the experimental dataEI on YMn 2 -based materials, where the transition is first- order. The 
present work amounts to a renormalization group derivation of some of the results obtained 
inil (staggered susceptibility), plus new results (specific heat coefficient, possible phase 
diagrams and the transition line equation) along with comparison to the experimental data 
for Ce-based metallic antiferromagnets undergoing second-order phase transition. It is also 
shown that generation of stiffness in the "soft" direction by short- distance fluctuations is 
negligible at all relevant momenta, and thus the quantum Lifshitz point studied in this paper 
is well-defined. 
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II. THE RESULTS 



The main results of this Section amount to establishing the leading low-temperature 
behavior of various quantities near a quantum Lifshitz point. This task is facilitated by the 
fact that the theory falls above its upper critical dimension. Moreover, since the stiffness in 
the "soft" direction, generated by the short-range fluctuations, turns out to be negligible (see 
Section VI), the thermodynamic and transport properties in the classical Gaussian region 
(T/v = 0) may be obtained from the Gaussian action with this stiffness set equal to zero: 

S{<P] = [ dr [ dx<\> a \S + \d T \ + Vjj + Vil (j> a , (1) 







where 6 is defined by the tuning parameter p and by the feedback of the quartic interaction 
in the Ginzburg-Landau action (||-fD as per 5 = (p — p c )/p c + const.T 5 / 4 . Thus the Neel 
temperature T N in this theory scales as T N ~ (p c — p) 4//5 . At p — p c , the correlation length 
£ii scales as T~ 5 / 8 , whereas the correlation length £j_ in the "soft" direction scales as T~ 5 / 16 . 

The leading exponent of the specific heat coefficient is given by the Gaussian contribution 
to the free energy F: 



an 



Q\\+Q± 



F = Tr log G 1 (q,uj) = J dz J d 2 q\\ J dqj_ coth — arct 
and yields the specific heat coefficient 

5C/T ~ T 1/4 . 

Resistivity due to scattering off anisotropic fluctuations can be estimated via the char- 
acteristic transport time Tt r given by the fluctuations of the magnetic order parameter near 
the transition. Since the critical fluctuations are antiferromagnetic and thus characterized 
by a finite wave-vector, the (1 — cos#) factor in the transport relaxation rate may be omitted 
in a qualitative estimate, which leads to 

q,u> 
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as realized indeed by the authors onllHla. Recently, it has been arguedEil that the observed 
resistivity p(T) ~ const. + y L2±ai may be a crossover phenomenon due to the interference 
between the impurity scattering and scattering by the critical fluctuations, conspiring to 
mimic a power-law behavior. At the moment it remains to be seen how sensitive this 
mechanism may be to the material-dependent factors such as the shape of the Fermi surface. 
At the same time, for a system with a true quantum Lifshitz point, the T 5//4 resistivity shall 
be accompanied by the C/T ~ T 1 / 4 scaling of the specific heat coefficient. 

In CePd2Si2 and CeNi2Ge2, the resistivity does exhibitliE a temperature exponent close 
to 5/4. However, the rest of the available data is less encouraging for the present theory. 
The Neel temperature scales linearly with pressure^ instead of obeying T/v ~ {p c — p) A ^. 
To my knowledge, the specific heat data on CePd2Si2 at the critical pressure is not yet 
available. The specific heat data on CeNi 2 Ge2 is still ambiguous, as both the C/T ~ T 1 / 2 
and C/T ~ logT have been reportedB 2 !. 

One can indeed find formal reasons (such as presence of the superconducting phase in 
CePd 2 Si 2 , or disorder in CeCu 5 . 9 Auo.i) why the Lifshitz point theory is not applicable to 
these materials. However, the T 1 ' 2 scaling of the resistivity in CePd2Si2 is seen at temper- 
atures up to 60 K, which is more than two orders of magnitude greater than the supercon- 
ducting transition temperature T c « 0.4 K. At the same time, virtually linear scaling of T/v 
with pressure persists up to about 5 K, which is again much greater than T c . These argu- 
ments (as well as recent finding0 of the u/T scaling and the inverse uniform susceptibility 
behaving as x 1 ~ const. + T°- 8±0A in CeCu 5 . 9 Auo.i) strongly suggest that the physics of 
these materials amounts to more than just a refined Gaussian theory above its upper critical 
dimension. 

III. THE SCALING EQUATIONS 

The assumption of an overdamped critical mode with the stiffness vanishing in one of 
the directions leads to the following effective action close to the quantum Lifshitz point: 
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S eff [4>]=S&>+S«> (2) 



o 



S (2) = dr dx<t> a 5 + \d T \ + Vjj + DVl + Vi Q (3) 



= [ dr [ dx \u{(i) a (i) a ) 2 + v^(V±(j> a ) 2 (j>p(j>p + v 2 ((f) a V \(p a ) 2 } . (4) 





with the frequency and the momentum cut-off set equal to unity. Here the mass term S and 
the stiffness D in the "soft" direction describe deviation from the quantum Lifshitz point, 
(3 = 1/T is the inverse temperature, Vn corresponds to the two "normal" directions, whereas 
Vj_ corresponds to the third, "soft" direction. The coupling constants Vi and v 2 describe 
the dispersion of the quartic coupling constant u\ they correspond to the only two possible 
linearly independent terms quadratic in Vj_ and quartic in <j) a , and describe generation of 
stiffness in the "soft" direction, as shown in Fig. [3]. 

The scaling variables in the theory are the mass term 5, the stiffness D in the "soft" 
direction, the temperature T and coupling constants u, V\ and v 2- To the lowest order in the 
latter three, the scaling equations can be derived e.g. by expanding the partition function Z 
to first order in then integrating out a thin shell near the cut-off, rescaling the variables 
and the fields, and then comparing the result with the original action^: 

dT(b) 



2T(b) (5) 

u{b) (6) 

viAb) (7) 
26(b) + 2{n + 2)ufx[T(b), 6(b), D{ 



dlnb 

du(b) 1 

dlnb " ~2 

dv 1>2 (b) _ 3 

dlnb ~ ~2 
d5(b) 
dlnb 

+ (nv 1 +v 2 )f 2 [T(b),6(b),D(b)} (8) 
^ = D(b) + (n Vl + v 2 )h[T(b),8(b),D(b)]. (9) 

Here n is the number of components of field <p a . The definitions of f\[T(b),8(b),D(b)] 
and f 2 [T(b),S(b), D(b)], as well as the details of the derivation, are given in the Appendix. 
Note that, following Millis@, I denote the running value of a scaling variable (e.g. T(b)) by 
indicating explicit dependence on the rescaling parameter b, whereas for the bare (physical) 



quantities the b dependence is omitted. Also note that truncating the equations at the first 
order in u, V\ and v 2 is equivalent to the assumption that the bare value of these couplings 
is small. 

An extra scaling variable D(b) in equations (0-0) compared with those of Millisi leads to 
appearance of two different ordered phases to be described in Section V. This sets the main 
distinction between a quantum Lifshitz point and a simple quantum critical point of Hertz 
and Millis, where there is only one ordered phase. 



IV. SOLUTION OF THE SCALING EQUATIONS 

In this Section, I obtain the qualitative solution of the scaling equations (§-|9|) following 
Millisi and noticing that /i[T(6), 5(b), D(b)} and f 2 [T(b), 5(b), D(b)] virtually do not depend 
on 5(b) or D(b) for 5(b), D(b) -C 1, and that both of them fall off rapidly as 5(b) exceeds 
unity. Thus one can neglect their dependence on 5(b) and D(b) for 5(b) < 1 and stop the 
scaling at 5(b) = 1. 

With these provisos, the formal solution of (H-0) reads 

T(b)=Tb 2 (10) 

u(b)=ub~ 1/2 (11) 

vi, 2 (b) = v lt2 b-V 2 (12) 

5(b) = 5b 2 + 2b 2 (n + 2)u [** ' dre~^ ' h[Te 2T ] 

Jo 

+ b 2 (n Vl + v 2 ) I n olTe-l r f 2 {Te 2T ] (13) 
Jo 

D(b) = b D + (nv x + v 2 ) r b dre-^fclTe 2 '' 

Jo 



(14) 



As in the case of a simple quantum critical point@, two major regimes exist, depending 
on whether the scaling stops when the running value of the temperature is much smaller or 
much greater than one. The first case is usually referred to as quantum while the second 
case is called classical. 
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To obtain the condition for the quantum regime, set T = on the right hand side of 
(|T3|), then integrate over r up to ln(b*) such that 5(b*) = 1, solve for b*, substitute b* into 
(|TCi|) and require T(b*) <C 1. The resulting condition is 



T < ri, n = -(n + 2)u/i[0] + ^(n«i + v 2 )f 2 [0]. 
5 7 



(15) 



If reversed, the inequality ( |15|) corresponds to the classical regime, where it is convenient 
to divide the scaling trajectory into two parts, corresponding to T(b) <C 1 (quantum) and 
Tib) 1 (classical). For T(b) 3> 1, one can estimate fi, 2 [T] as fi[T] ~ BiT, f2[T] ~ -B 2 T. 
Then, for T(6) ^> 1, the equations (^-§) may be recast in terms of new variables U(b) = 
u{b)T{b),V 1>2 {b)=v 1;2 (b)T(b): 

dT(b) 



dlnb 

dU(b) _ 3 
dlnb ~ 2 
d^i,2 (6) = 1 
dlnfc 2 
dS(b) 

dlnb 
dD(b) 

dlnb 



2T{b) 
U(b) 

25(b) + 2B 1 (n + 2)U(b) + 5 2 (nl4(6) + V 2 {b)) 
D{b) + B 1 {nV 1 {b) +V 2 {b)). 



(16) 
(17) 
(18) 
(19) 
(20) 



The initial conditions correspond to b = b such that T(b) — 1 and read 

T(5) = 1 

U(b) = u(b) = uT 1/4 
V 1>2 {b) = v li2 T 3 / A 

5(b) = - fn + A(n + 2)uT 5/4 



D(5) 



r 2 + -(nvi + v 2 )T 5/4 



(21) 
(22) 
(23) 
(24) 
(25) 



where r 2 is defined by 



r 2 = J D + -(nu 1 + u 2 )/ 1 [0]. 



(26) 
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For T 3> ri, both £7(6) and Vi^fe) are small, which justifies using linearized equations near 
b — b. Neglecting the higher powers of T and assuming that u and Dip are of the same order 
of magnitude, the solution of (|T6| - pOD is 

T{b) = T{b)b 2 (27) 

U(b) = U(b)b 3/2 (28) 

V h2 (b) = V 1>2 {b)b 1 / 2 (29) 

5(b) = b 2 [5(b) + AB lU {n + 2)T 1/4 ] (30) 

D(b) = b [D(b) + 2B 1 (nv 1 + v 2 )T 3/ *] . (31) 

Now one can ensure the consistency of using the linearized equations by demanding that, 
when scaling stops at S(b*) = 1, the coupling constant U is small: 

U(b*) < 1. (32) 

This condition corresponds to the Ginzburg criterion^ and is violated only very close to the 
line S(b*) = 0, which defines the Neel temperature T/v as a function of u and r\: 

n + (A + AB X ) (n + 2)uT b J A = 0. (33) 

The Ginzburg criterion fl3l|) is violated in a narrow window 5Tn/Tn ~ u 

1/3 T V12 <<; x of 

strong classical fluctuations. 

To establish connection with experiment, it is assumed that T\ and r 2 are both propor- 
tional to the deviation (p — p c ) of the control parameter p from its critical value p c . This is a 
reasonable assumption given that the theory at hand is above the upper critical dimension. 
Thus T/v scales as per 

T N ~(p c -p) 4 / 5 . (34) 
However, which two phases does the line pi) separate? 
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FIG. 1. Schematic phase diagram in the vicinity of a Lifshitz point. The ordered phase I 
corresponds to ordering at the wave vector q = and is shaded by horizontal lines, whereas the 



ordered phase II corresponds to ordering at q = ±^—D/2 at D < and is marked by zig-zag 
shading. 

V. THE PHASE DIAGRAM 

To answer this question, one has to recall that, as mentioned above, Lifshitz point is a 
tricritical point. Thus, generally, one shall expect to see two different ordered phases and 
one disordered phase in its vicinity, as illustrated by the following toy-model expression for 
the free energy F at a finite-temperature Lifshitz point: 

F ~ [6 + Dq 2 + g 4 ] 0. (35) 

As follows from (|35|), the region (D > 0, 5 < 0) corresponds to the phase ordered at q — 0, 



whereas in the region (D < 0, 5 < (D/2) 2 ) one finds ordering at wave vectors q = ±y—D/2. 
This phase diagram is shown in Fig. |T]. 

To establish the phase diagram near a quantum Lifshitz point, one shall draw the curves 
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FIG. 2. Schematic phase diagram in the vicinity of a quantum Lifshitz point. As in Fig. |], 
horizontal shading denotes "commensurate" ordering (at q = 0), whereas zig-zag shading corre- 
sponds to "incommensurate" order (at q = iy— D/2). Wherever shown, dashed line corresponds 
to 5 = (D/2) 2 at D < 0. The bold lines defined by D = and 5 = are determined by equations 
© and @. 

of 5(6*) = and D(b*) = in the (T — p) plane and compare with Fig. [p. Mutual position 
of the two curves depends on the relative magnitude of u, V\ and t>2, and on the signs of 
the proportionality coefficients between r 12 and {p — p c )- Six possibilities arise, as shown 
on Fig. ^. As in Fig. [I], horizontal shading denotes "commensurate" ordering (at q — 0), 
whereas zig-zag shading corresponds to "incommensurate" order (at q = ±J —D/2). Case 
(d) can be ruled out on physical grounds, as it corresponds to ordering at high temperatures. 
The possibilities illustrated in Fig. (a) and (b) are irrelevant to the results obtained for 
CePd2Si2 or for CeCue-xAu^, as they lead to existence of two ordered states, whereas no 
trace of second transition has been found in either of the materials of interest. The remaining 
cases (c), (e) and (f) are the most interesting for us, as they show how a quantum Lifshitz 
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point can mimic a "regular" quantum critical point with only one ordered phase. Case 
(f) corresponds to "commensurate" ordering, while (c) and (e) describe "incommensurate" 
order. In the latter two cases, the transition line between the "incommensurate" and the 
disordered phases is described by the equation 5(b) = (D(b)/2) 2 , shown in Fig. Q by dashed 
line, with 5(b) and D(b) given by ( |5U| - pi~D . At low temperatures, this line asymptotically 



coincides with the line 5(b) = given by fl34j) . 

Currently, the structure of magnetic order in CeC^-zAu^ is being mapped out exper- 
imentally by several groups@. However, if a Lifshitz point described above is realized in 
CeCu6-a;Au x , the exact character of ordering is irrelevant for the physical properties in the 
"classical Gaussian" region, roughly corresponding to the vicinity of the line (T > 0, T c = 0). 

VI. GENERATION OF STIFFNESS 

As shown in Section III (see Fig. |3|), near the quantum Lifshitz point the short-range 
fluctuations do generate stiffness in the "soft" direction, even though it is absent exactly at 
the critical point. To check the importance of this effect at finite temperature, one has to 
compare the generated quadratic term Dq\ with the quartic term q\. The comparison can 
be done easily using the solution (0-0) °f the scaling equations 

At each step of the renormalization procedure, one focuses on momenta and frequencies 
near the current value of the cut-off. Recalling the agreement to set the cut-off equal to 
unity, one finds that rescaling factor b corresponds to frequency and momenta q\\(b) = 1/b, 
q± (b) = 1/Vbandu(b) = 1/b 2 . 

Thus the relative importance of the generated term D(b)q\(b) is given by its comparison 
with the mass term 5(b) and with the quartic term q\(b). To do the comparison in the 
classical Gaussian region which we are studying, one shall set r± = r2 = and then, using 
( [llj - Il ), estimate 5(b) and D(b) at a running value of b. 



In the regime of "quantum" renormalization (T(b) <C 1), the sought estimate is given by 

5(b) ~ A(n + 2)ub~ 1/2 , 
14 




Q = A(p,q);\J=B(p) 



STIFFNESS: CT) <=l ' CD 



FIG. 3. Generation of stiffness in the "soft" direction by the renormalization process. Prop- 
agators of "fast" degrees of freedom, integrated out at each step of renormalization, are shown 
by dashed lines. The free energy correction in the second order in point-like quartic interaction u 
(black dot) generates dispersion of the effective quartic interaction (shaded circle). The latter, in 
turn, gives rise to the stiffness in the "soft" direction (shaded square). 

D(b) ~ A(n Vl +v 2 )b^ 2 , 
D{b)q 2 ± {b) ~ A{n Vl + v 2 )b-^ 2 . 

qi(b) ~gj(&) ~ M(6) ~r 2 . 

Applicability of our scaling equations requires weak coupling (u,v i; 2 <C 1). Thus, D{b)q 2 ± {b) 
is negligible compared with q\(b) at any value of b up to the point when 5{b) and q\{b) 
become of the same order of magnitude. Therefore, the stiffness generated by short-distance 
fluctuations under the renormalization flow can be safely neglected in our study, which in 
turn means that the quantum Lifshitz point is well-defined. 
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VII. CONCLUSIONS 



In this paper, I studied the classical Gaussian region near a quantum Lifshitz point. 
This corresponds to defining the low-temperature behavior at the point where the Neel 
temperature and the stiffness in the "soft" direction simultaneously become equal to zero. 
The Neel temperature was found to scale as T/v ~ (p c — p) 4//5 . The specific heat coefficient 
was found to have a T 1 / 4 anomaly, whereas the resistivity was shown to exhibit T 5 / 4 scaling. 
Of these results, only the resistivity exponent finds experimental support (in CePd2Si2 and 
CeNi 2 Ge2), while the predictions for the shape of the transition line appears to fail in fitting 
the experimental data. The situation with the specific heat data is not yet entirely certain. 

Regardless of possible reasons, discussed briefly in Section II, the failure of simple phe- 
nomenology based on the assumption of a quantum Lifshitz point appears to be unambiguous 
for CePd 2 Si 2 and CeCug-xAuj;. More generally, it appears that the correct theory of the 
transition shall fall below the upper critical dimension, while the Hertz-Millis theory and its 
present refinement are all essentially Gaussian. 

At the moment, the source of the puzzling behavior of these Ce-based materials near a 
zero-temperature transition remains unclear. A comprehensive experimental study (specific 
heat, neutron scattering, thermal transport and NMR/NQR/yuSR) might help to resolve 
some of the pressing issues. An intriguing possible direction of theoretical research would 
be to study the yet poorly understood interplay between the incipient correlations in the 
conduction sea and the development of the Kondo effectll. 

I am indebted to G. Aeppli, P. Coleman, A. Millis, D. Morr, A. Rosch and A. Schroder 
for discussions related to this article, and to Y. Aoki and T. Fukuhara for discussions of their 
data. The work was started at Rutgers University, where it was supported by the National 
Science Foundation (grant number DMR-96-14999). Work at the University of Illinois was 
supported in part by the MacArthur Chair endowed by the John D. and Catherine T. 
MacArthur Foundation at the University of Illinois. 



16 



APPENDIX 



In this Appendix, I derive the scaling equations ([5|-§). The first term on the right hand 
side of each equation corresponds to rescaling of the variable under an infinitesimal time- 
and length-scale transformation. It can be obtained e.g. by following HertJ. Rewrite the 
action (|§-fD i n the momentum and frequency domain, replacing the sum over the Matsubara 
frequencies by an integral up to the cut-off (set equal to 1): 



S ef f[<P}=S^ + S^ 

s{2) r^du f 1 dql 
S (4) 



2ir J (2vr) s 



dq± 
2vr 



5 +\uo\ + q\ + Dq\ + q\ 



(36) 
(37) 



1 duj f 1 dq 2 ri dq± 
2^ J (2tt) 2 J 



-i 3 



J 1.2,3 



U{<p a <p a f + M<l±4>a) 2 (p(3(pf3 + V 2 {<l> a q±.<j> a )' 



(38) 



Then integrate out a thin shell between the original cut-off (equal to 1) and the new 
cut-off (equal to 1/6) in the q\\ space. Now define thin shells in the ui and q±_ spaces in 
such a way that upon being integrated out, they would admit rescaling of all the variables 
(u,qn,q±,S,D,u,Vi,v 2 and <p a ) in such a way as to bring the remaining part of back 
exactly to the form (18) with the new values of 5, D, u, V\ and v 2 - The only choice of rescaling 
factors which allows this corresponds to rescaling qn by 1/b, integrating out u between 1 
and 1/b 2 and rescaling it by l/b 2 , integrating out q±_ between 1 and 1/y/b and rescaling it 
by 1/y/b - at the expense of rescaling 5 and T by b 2 , D by b, u by 1/y/b and t>i >2 by 6~ 3//2 , 
which corresponds precisely to the first terms on the right hand side of (H-^|). 

Now, I will obtain the remaining terms in (|-|9|). First, expand the partition function Z 
to first order in S^: 



Z () 



(S (4) > 



Using Wick's theorem^, the average (S^) can be conveniently rewritten as 

(S^) = n(n + 2)u(<P 2 ) 2 + n{nv x + v 2 ) (0 2 > (q 2 ^ 2 ) , 



where ((f) ) and {q±(p ) are defined as per 
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l2\ 



1 dz ( z \ f 1 dqf\ f 1 dq 



7T 



Z 



) = L - c ° th b^Tn 77^ "^ -5 I rim I .2 I F^2 I ^ (39) 



2T(6) y/ 7 (2vr) 2 7 2vr z 2 + [5(6) + gj + D(6)g 2 + <ffi< 



v coth (^)J^ ^r ^ + [^) +9 g + D(%i + gip - (40) 

The next step amounts to integrating out a thin shell near the cut-off in the expression for 
(S^). As described above, such a shell has width 1 — 1/6 in the q\\ direction, 1 — 1/6 2 in the 
z direction and 1 — 1/ Vb in the q± direction. This integration generates the sought terms 
in (BE), with 



1 dq L 



2T(b) J J 2?r z 2 + [5(b) + 1 + D(b)q\ + 



2 ,J 1 \ f 1 d 4 f 1 dq. 
H — cotn 



2T(b) J J (2tt) 2 J 2tt 1 + [5(b) +qf + D(b)q\ + q\f 
1 /-i dz f z \ r 1 dqj z 
Wo ~ COt \mp)J J (27r) 2 z 2 + [5(6) + g 2 + J D(6) + l] 2 ' 

(41) 



2tt z* + [5(6) + 1 + D(b)q\ + gj] 2 



+ 2 coth / i ^ r H i 



71 



2T(b) J J (2tt) 2 J 2vr 1 + [5(b) + q\ + D(6)g 2 + g*] 2 



1 f 1 dz ( z \ r 1 d i z i 



H / — coth 



4tt7o tt \2T(b) J J (27r) 2 z 2 + [5(6) + g 2 + D(6) + l] 2 ' 

(42) 



thus completing the derivation of the scaling equations (|-|9|) 
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